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1. INTRODUCTION 

The economic load dispatch (ELD) of generation is one of the most crucial tasks in modern power 
systems. It promises better dispatch scheduling to mitigate the increasing cost of fuel for various types of 
thermal power plant. Solving the ELD problem is a potentially significant economic solution to power system 
planning and operations. ELD’s goal is to commit the required power consumption by planning the different 
types of power generation units such that secure total cost production is the cheapest possible while fulfilling 
the system’s equality and inequality constraints. It keeps the cost of producing electricity to a minimum price 
by properly allocating workloads among the generating units of the plants with various operating efficiencies, 
type of fuel cost and total transmission losses of systems. This optimum solution to the problem of generating 
power contributes significant economic benefits to the operation of the power plant. 

Initially, the ELD problem was solved using traditional techniques such as linear, quadratic and 
nonlinear programming methods as provided in [1-5]. The conventional techniques have a higher 
probability of becoming trapped in local minima solution due to the complex fuel cost function problem 
related to the highly nonlinear characteristic of present power generating units such as ramp rate limit and 
prohibited operating zones. Nowadays, in order to handle the nonlinear fuel cost function, many advanced 
optimization techniques based on nature-inspired meta-heuristic has been implemented on ELD such as 
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Genetic Algorithm [6], Particle Swarm Optimization [7], Artificial Bee Colony [8, 9], and Ant Colony 
Optimization [10], Cuckoo Search Algorithm [11, 12], to name a few. 

In 2014, Mohseni, Gholami, Zarei and Zadeh [13] introduced a new meta-heuristic algorithm called 
Competition Over Resources (COR), inspired by a group of animal communities which compete for 
resources. The COR algorithm was evaluated for quality results by addressing various literature optimization 
issues in [14, 15]. From the literature, the COR algorithm is used to evaluate eight different benchmark 
features and concludes that this optimization method aims to be a global solution better than the PSO and GA 
algorithms. COR also has been tested for its capacity to address the ELD problem in [16] to find the optimal 
generating costs on 6-unit and 15-unit test systems. 

One crucial problem while applying Optimization to the ELD problem is that the solutions will not 
always satisfy the inequality and equality constraints simultaneously. Most ELD constraint optimization 
problems have adopted the Penalty Function Strategy (PFS) approach [17-21] to handling constraints because 
of its simple implementation. The penalty function approach involves several penalty parameters which need 
carefully determined tuning value to obtain a feasible solution. This approach needs extensive 
experimentation for setting up the appropriate parameters needed to define the correct penalty function 
parameter. This paper proposed a Feasible Solution Strategy (FSS), which prevents the creation of the 
infeasible solution. Using the FSS method, the infeasible particle will be relocated to new particles until a 
feasible solution condition is satisfied. 

This paper compares the performance of COR using PFS and FSS to achieve a quality solution for 
the ELD problem. Due to simulation based on practical operating conditions of the power system, 
the equality and non-equality constraints is considered, which increased complexity in search for a quality 
solution. The objective of using PFS and FSS with COR optimization to overcome the ELD issue is to 
compare the performance with other optimization techniques in terms of its quality as a global solution and 
convergence capability. 


2. ELD PROBLEM FORMULATION 

The ELD issue, which deals extensively with financial concerns, refers to the amount of power to be 
assigned from all generators in an attempt to minimize the cost of supplying the required electrical energy 
subject to several constraints of multiple generating units. The ELD’s mathematical optimization has three 
primary parts, which are the optimization variable of the problem, the goal of the objective function and 
constraints. 


2.1. Optimization Variable 
The optimized input variable is the real output power of generating units, 


P = [P4,Po,...,Pwg| (1) 
Where, Ng is the total number of generating unit in the power system. 
2.2. Objective Function 

The objective function of ELD is to supply power via each generator unit for the request demand 


load with minimum generation total fuel cost. ELD’s objective function can be expressed through the 
problem of mathematical optimization, 


minimize F., = SF, (F ) (2) 
j=l 
Where 
F;(P;) = generation cost of the j" generator ($/h) 
P; = power of the j” generator (MW) 


Moreover, the cost coefficient of each generator stated as, 


F,(P,)=a,+B,P,+7,P; (3) 


J 


Where, 
a;, 8; and y; = coefficients of the 7" generating unit 
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2.3. Constraints 
To achieve ELD's optimum value, the objective function has been subject to the following limitation 
constraints: 


2.3.1 Power Demand Constraint 
Total power generation must be satisfied with the total power demand and power losses as follows, 


Ng 4 
DP =P, +P, (4) 
j=l 
Where, 
Ng 
pie = Total power generation (MW) 
J 
jel 
Pp = Total power demand (MW) 
Pr = Total power transmission loss (MW) using Kron’s loss formula 
The total power transmission losses calculated as follows, 
Ng Ng Ng (5) 
P= SD Pae + DP,Bio + Boo 
jal k=l jel 
Where, 
Bj» Bjo, Boo = B-coefficients or loss coefficients. 
2.3.2 Generation Constraint 
Each generator's power generation must be within its operating limits as follows, 
| ial a oe eo (6) 


Where, P}””" and P;"“ is minimum and maximum generation limit in MW. 


2.3.3 Prohibited Operating Zones 
The system involves certain restricted areas related to physical component constraints as follows, 


min I 
Proz.j = Proz.j = Pro. 
uu 1 =, 
Prot, j.t-1 = Proz,j = Proz, ja» k ~ 2, oe : Mj 
u max 
Proz.ian = Proz.j Ss Proz.j (7) 


Where, 

Prozjk = minimum border of k’" restricted zone of the j” generator 
Proz, ike = maximum border of k" restricted zone of the j" generator 
nj = number of restricted zones of the j” generator 


2.3.4 Ramp Rate Limits 
The operating range of all operating units is restricted by their ramp rate limits in the actual power 
generation process for each unit. An increase or decrease of power generation is as follows, 


Prert.j = Pak s UR; 
Pigs = Pert, j s DR, 
max Pig Pye? SDR) < Pigg Sin Bay Pi? OR; ) (8) 
Where, 
UR; = up ramp limit of the j” generator (MW/h) 
DR; = down ramp limit of the j” generator (MW/h) 
Prrij’ = previous output power of the j” generator (MW/h) 
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3. COR OPTIMIZATION ALGORITHM FOR ELD 

The COR begins at each year with the groups divided based on the food supply for each territory. 
At this stage, each group is looking for food in its territory. An active participant will periodically remove the 
weaker participant from the same region throughout the process so that the final process will consist of 
excellent competitors only. When food sources in a particular field are reduced, the members of the group 
will also be decreased. Members possibly will migrate to join a group in regions where there is plenty of food 
resources. Finally, the group with fewer resources is eliminated at the end of the year, leaving only the group 
with the most food sources. Figure 1 shows the pseudo-code of the COR optimization algorithm 
implementation to solve the ELD problem. 





The Initialisation of COR parameters 
Evaluation of generation costs with constraint handling technique 
Select each group's best agent 
while (Condition of termination is not fulfilled) do 
Evaluate the distance among group agent 
Generate random search inner and outer agent 
Evaluate new generation cost with constraint handling technique 
Update each groups' best agent 
Increase an agent from the greatest group and eliminate an agent from the weakest group 
if (group population = d,a:c) 
Remove the weakest group and split the greatest group into two new group 
end if 
end while 











Figure 1. Pseudocode for COR optimization algorithm 


3.1. Initialisation 

Identify parameters of COR such as minimum (P,,;,) and maximum (Prax) generation limit of each 
unit, number of iterations (Ni;e-), number of population (N,.,), number of groups (Neroup), rate of death (dae) 
outer searching factor (D,.,) and population ratio between inner/outer neighborhood (P,,;,). The generator’s 
active power generation is identified as an ELD problem input variable. Then, the population was randomly 
distributed, ranging over their maximum and minimum limit according to the generator limit stated in (6). 


3.2. Evaluation of Generation Cost with Constraint Handling Technique 

The fitness of each population is assessed with the objective function in (2) using the Penalty Factor 
Strategy (PFS) or Feasible Solution Strategy (FSS). The result of the fitness cost calculation is graded from 
the highest to the lowest solution. Then, the agents were equally distributed into the number of groups, Neroup 
and the best solution was identified as the best group agent for each group. 


3.2.1 Penalty Factor Strategy (PFS) 

This strategy penalizes ineffective solutions by multiplying a persistent penalty for these alternatives 
that violate the limitations. To satisfy constraint limitation and disallow the workable area, the individual 
population’s fitness performance is measured using (9), which is the combination of objective function 
equation between (2) and PFS linked with related constraints. This evaluation feature is used to obtain the 
smallest cost generation value while fulfilling the limitation of equality and non-equality constraint problem, 
as stated in (4-8). 


Ng 


FR)= DFP) +k 


Ng 


DP. - (+h) 


j=l 


Ng Ng (9) 
+k, > Pio yr k, Pee: j 


jal joi 








Where: 
k1, k2, k3 = constraint constant value 


3.2.2 Feasible Solution Strategy (FSS) 

The PFS technique does not ensure that all solutions generated by the optimization algorithm meet 
the limitations of equality constraint because it is exceedingly difficult to obtain when considering nonlinear 
problems in ELD. Using FSS, handling constraints based on the repair of the unfeasible solution is 
implemented to guarantee that all the solutions generated are satisfactory through the optimization process. 
The details of FSS procedure is shown in Figure 2. 
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Step 1: calculate power balance error (4P) using 
AP = Pp - (X(Pc) + Pi) 
where: 
Pp = power demand 
Pg = power generation 
P,  =transmission losses 
Step 2: randomly select generator using 
n = fix((rand*d) + 1) 
where d is the number of generators 
Step 3: while (|4P| > 0.00001) 
Updated new generation for select generator using 
Po(n) = PG(n) - |AP| 
Check P,(n) if exceeds the maximum or minimum generator operating limits 
then assign value to its limits 
Check P,(n) if located between trapped area POZ then assign value to nearest 
it is boundaries 
Calculate 4P using the new value of Pg and P;, 
Randomly select generator n 
end while 
Step 4: Update value of Pg 











Figure 2. Constraint handling based on FSS 


3.3. Determination of Territory 

Each group’s territory is described by using the Euclidean distance among the most excellent agents 
in the group. The minimum value of the distance between territories, d”” is the option rate to specify inner 
territory. The maximum and minimum power generation boundary for inner territory determined in (10), 


P. 


inner 


P. 


inner 


min =P. —~q™ 
max a +qmu (10) 


The maximum and minimum power generation for outer territory determine in (11), 


™ = P,— (dy!) 
max = P.+(d,,,x1) (11) 


outer 


outer 


Where / is the different value between Pig, and Pynin and ds-, is the option rate of the outer space search 
capability between 0 and 1. 


3.4. New Population Generation and Evaluation 

The population quantity between the inner and outer territories of each group must be determined 
using the ratio of inner and outer neighbourhoods, P,.,. Some quantity of the agent’s new population in the 
group was produced randomly using the inner territory to locate a prospective optimum area within the group 
boundary. The left agent is obtained randomly using the identified outer territory with a larger searching 
space area to raise the opportunities of a random agent discovering optimal region outside group margins. 

Finally, the new generated agents are set to the maximum and minimum boundary shown in (6). 
Then, all agents are evaluated using the process details in step 3.2 and update the position of the new best 
group agent. 


3.5. Update Group Members 

The group with the highest results from one of the group members will add a new member to its 
community while the poorest-performing group will eliminate one of its members. The competitive 
organization will increase its population from this process, to discover potential wealthy assets and the least 
productive organization will reduce its members. 
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3.6. Termination Condition 

This ELD optimization problem uses the number of iterations, Nj, as an ending criterion for the 
process. The COR algorithm ends the process when the highest iteration is accomplished, and the optimal 
result is recorded. Else, this process in step 3.3 is repeated. 


4. RESULTS AND ANALYSIS 

The PFS and FSS constraints handling approaches for the ELD problem have been tested using the 
COR algorithm on two standard different test systems in order to compare their performances. All simulation 
work is carried out with MATLAB programming on a PC specification with an Intel Core 17 Dual-core 
2.80GHz, 8GB RAM. For all test cases, the population size is set to 100, 200 maximum iterations and run 30 
times, respectively. The parameters setting for COR used for both case studies are listed in Table 1. 


4.1. 6-Unit Test System 

This case study comprises six generators together with constraints of ramp rate limits, real power 
balanced, prohibited operating zones, and generating limits. Total active power demand is 1263MW, 
and transmission loss is accounted in this test. The parameter of this test system is obtained from [22], 
and the results are compared for the cases of GA [22], PSO [22], and NPSO-LRS [23]. 

Table | demonstrates the statistical results attained after 30 runs of PFS and FSS. The result shows 
that the COR-FSS method can produce a better solution on the minimum cost compared with COR-PFS and 
other optimization techniques. Except for GA, all the other algorithms give almost the same minimum 
generation cost of COR-FSS. Both COR-PFS and COR-FSS satisfied the system constraints. 

For robustness test, the dominance of the COR-FSS over COR-PFS, GA and PSO can be observed. 
The maximum cost and average cost values obtained by COR-FSS are very close to the minimum cost value 
and the lowest value of the standard deviation (SD). Result proves that COR-FSS is more robust and 
constitutes a more consistent solution compared to other techniques. 


Table 1. Optimal Results for the 6-Unit Test System 
Optimization Method 





Power Generation (MW) 





GA PSO NPSO-LRS COR-PFS COR-FSS 
G, 474.81 447.50 446.96 448.0018 447.4870 
G 178.64 173.32 173.39 173.0716 173.3266 
G; 262.21 263.47 262.34 264.6342 263.4636 
Gy 134.28 139.06 139.51 138.4441 139.0714 
Gs 151.90 165.47 164.71 165.4049 165.4668 
Go 74.18 87.13 87.98 86.4171 87.1427 
Total Power Generated 1,276.03 1,276.01 1,275.94 1,275.9737 1,275.9586 
Power Losses (Ploss) 13.02 12.96 12.94 12.9737 12.9586 
Total Generation Cost ($/h) 15,459 15,450 15,450 15,449.9235 15,449.8994 
Minimum Cost 15,459 15,450 - 15,449.9235 15,449.8994 
Maximum Cost 15,524 15,492 - 15,478.9031 15,449.8995 
Average Cost 15,469 15,454 - 15,460.6548 15,449.8994 
SD 0.0570 0.0002 - 7.5517 3.3043e-05 





Figure 3 shows the convergence performance between COR-FSS and COR-PFS for achieving the 
best solution cost. It shows that the COR-FSS strategy reaches the lowest cost of generation quicker than 
COR-PFS. COR-FSS generates the lowest cost in the initial iteration since it can eliminate the conditional of 
constraint generation compared to COR-PFS, which is affected by the penalty factor for not exceeding the 
conditional of constraint. 
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Figure 3. Convergence behaviour of COR-FSS and COR-PFS for 6-unit test systems 


4.2. 15-Unit Test System 

This case study consists of 15 generating units of large-scale test systems with a total active load 
demand of 2630MW. It also considers the ramp rate limits, real power balanced with transmission losses, 
prohibited operating zones, and generating limits. The parameters of this test system are obtained from [22], 
and the results are compared with the cases of GA [22], PSO [22], ABC [24] and GA-API [25]. 

The statistical results of maximum, minimum and average cost achieved after 30 individual runs of 
PFS and FSS are based on the COR algorithm and compared with other methods as shown in Table 2. 
The result shows that COR-FSS obtained the best global solution with the value of 32,704.4499 $/h, followed 
by COR-PFS by 32,717.0105 $/h. Both techniques satisfied the limitation set by the test system constraints. 
It observed that the COR-FSS technique could produce better quality solution generation cost and the lowest 
value of standard deviation among all stated technique. 

The COR-FSS technique produced the lowest standard deviation (SD) with the value of 3.3154x10* 
compared to COR-PFS. However, COR-PFS’ value of standard deviation is higher than GA, PSO, ABC and 
GA-API. COR-PFS demonstrated its capability to obtain robust and consistent minimum results of cost 
generation solution. 

Figure 4 shows the convergence behaviour of COR-PFS and COR-FSS for a maximum iteration of 
200. COR-FSS achieved the lowest cost generation in the early iteration of 20 and gradually move towards 
optimal results by the end of iteration, which is much faster compared with COR-PFS. 


Table 2. Optimal Results for the 15-Unit Test System 
Optimization Method 








Pane Geena) GA PSO ABC GA-API__ COR-PFS__COR-FSS 

G, 415.31 439.11 454.2778 454.70 454.7253 455.0000 

Go 359.72 407.97 369.7131 380.00 379.3048 380.0000 

G; 104.43 119.63 130.0000 130.00 129.9789 130.0000 

G4 74.99 129.99 124.3210 129.53 129.7756 130.0000 

Gs 380.28 151.07 163.1341 170.00 169.1184 170.0000 

Ge 426.79 460.00 460.0000 460.00 459.9703 460.0000 

G, 341.32 425.56 405.4317 429.71 430.0000 430.0000 

Gg 124.79 98.57 85.6483 73.35 103.7853 71.7561 

Go 133.14 113.49 92.1289 34.96 48.9414 58.9054 

Gio 89.26 101.11 157.4626 160.00 133.7865 160.0000 

Gu 60.06 33.91 74.5293 79.75 79.2430 80.0000 

Gi 50.00 79.96 79.8057 80.00 79.8047 80.0000 

Gi3 38.77 25.00 25.0000 34.21 25.0048 25.0000 

Gis 41.94 41.41 19.3117 21.14 18.2376 15.0000 

Gis 22.64 35.61 20.8153 21.02 18.7265 15.0000 

Total Power Generated 2,668.40 2,662.40 2,661.5795 2,660.36 2,660.4030 2,660.6615 
Power Losses (Ploss) 38.28 32.43 31.5795 30.36 30.4042 30.6615 
Total Generation Cost ($/h) 33,113 32,858 32,787.8365 32,732.95 32,717.010 32,704.4499 
Minimum Cost 33,113 32,858 32,787.8365 - 32,717.0105 32,704.4499 
Maximum Cost 33,337 33,331 - - 32,981.0896 2,704.4535 
Average Cost 33,228 33,039 32,791.5366 - 32,794.4522 32,704.4503 
SD 0.0087 0.0070 2.4746 - 44.2100 3.3154e-4 
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Figure 4. Convergence behaviour of COR-FSS and COR-PFS for the 15-unit test systems 


5. CONCLUSION 

This paper comparatively studies the performance of two types of constraint handling techniques, 
which are PFS and FSS on COR optimization technique to solve the issue of non-convex ELD problem 
considering generator limit, power balanced, ramp rate limits and prohibited zones. The two test system is 
used to testify the effectiveness of the proposed method. Based on this research, the constraint handling 
technique affects the optimization algorithm efficiency of finding better quality solutions. Comparing the 
outcomes acquired from the COR-PFS and COR-FSS demonstrated that the COR-FSS technique is highly 
efficient in continuously offering superior solutions in term of searching for a global solution, 
quick convergence and robustness for non-convex ELD issues compared to COR-PFS and other optimization 
techniques. 
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